Differential-thermal analysis around and below the critical temperature Tc of various 

low-Tc superconductors: A comparative study 
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We present specific-heat data on the type-II superconductors VsSi, LuNi2B2C and NbSe2 which 
were acquired with a low-temperature thermal analysis (DTA) technique. We compare our data 
with available literature data on these superconductors. In the first part we show that the DTA 
technique allows for fast measurements while providing a very high resolution on the temperature 
scale. Sharp features in the specific heat such as at the one at the transition to superconductivity are 
resolved virtually without instrumental broadening. In the second part we investigate the magnetic- 
field dependence of the specific heats of VaSi and LuNi2B2C at a fixed temperature T — 7.5 K 
to demonstrate that DTA techniques also allow for sufficiently precise measurements of absolute 
values of Cp even in the absence of a sharp phase transition. The corresponding data for VsSi and 
LuNi2B2C are briefly discussed. 

PACS numbers: 65.40.Ba, 74.25. Bt, 74.25.Dw, 74.70.Ad, 74.70.Dd, 75.40.Cx 



INTRODUCTION 



The specific heat is a bulk thermodynamic quantity 
determined uniquely for every material by its spectrum 
of excitations. The measurement of the specific heat is 
a basic technique to reveal the physical properties of a 
material because it can, in principle, be calculated ab 
initio from a suitable physical model. The prediction 
of a linear specific heat, for example, is one of the most 
important consequences of Fermi-Dirac statistics for 
electrons in a metal, and its measurement provides a 
simple test to the electron gas theory of a metal. Because 
the specific heat is a quantity that is probing the whole 
sample volume, it is significant to characterize volume 
effects such as the superconducting state in a material. 
Therefore, specific-heat measurements have traditionally 
been of great importance for investigations in the field 
of superconductivity. The electronic specific heat of 
superconductors can be expected to yield information 
about the nature of the superconducting state; and its 
temperature (T) dependence should in particular be 
related to the energy gap [1]. Specific-heat measure- 
ments can also provide information about fluctuation 
effects [2-6], and they are particularly suited to measure 
phase transitions and to explore the phase diagram of 
superconductors. Improvements in the techniques for 
specific-heat measurements are therefore of special value. 
As an example, the first-order nature of the melting of 
the flux-line lattice in high-temperature superconductors 
in the mixed state has been unambiguously proven for 
the first time by high-resolution specific-heat measure- 
ments based on a differential-thermal analysis (DTA) 
method developed for this purpose [7, 8]. The DTA 



method used in our laboratory is particularly suited to 
detect sharp phase transitions, e.g., first-order phase 
transitions, and it does not require large sample masses. 
Since high-quality samples often only exist in form of 
small crystals, such a sensitive method is of particular 
interest. In addition, it is often desirable to choose a 
method that yields a high data-point density within 
a reasonable measuring time. Common techniques, 
such as standard hcat-pulsc or relaxation methods, are 
sensitive but very time consuming. As an example, 
the relaxation method implemented in a commercial 
PPMS platform (Quantum Design) usually takes several 
minutes per data point for data acquisition, while the 
DTA technique takes only approximately 3 s per data 
point, at an impressive data-point density of typically 
170 data points per Kelvin. High resolution and high 
data-point density are needed for the observation of 
small and sharp effects. In the first part of this work 
we are comparing corresponding DTA measurements on 
VsSi, LuNi2B2C, and NbSe2 with measurements done 
in a commercial PPMS platform and with data from the 
literature, with a special focus on the transition width 
of the normal to the superconducting state. 

Historically, specific- heat (cp) measurements using 
DTA techniques have been proven to be particularly 
powerful when studying relative changes in Cp at sharp 
phase transitions. However, to extract exact absolute 
values, or to identify smooth, featureless trends (such as 
the exact functional temperature dependence of Cp), is 
not easy using this technique due to various reasons [7]. 
In superconductors, not only the discontinuity in the 
specific heat at Tc contains valuable information about 
the superconducting state, but also the temperature 
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and the field dependence of this quantity below Tc are 
of utmost interest. From the temperature dependence 
Cp(T) one can draw conclusions about zeroes in the 
gap function, which may help clarify questions about 
the nature of the order parameter [9]. A power-law 
dependence of Cp on temperature has been observed, 
for example, in YBa2Cu307 [10], with a term for a 
magnetic field H = and an H^/^T term for iJ 7^ 
and low T, with a crossover to a stronger T dependence 
at high T. Peculiar T dependencies, very distinct from 
the standard one-band ,s-wave exponential law, can also 
be expected in multigap superconductors, as has been 
measured for example in MgB2 [11]. On the other hand, 
the field dependence of Cp{H) at constant temperature 
may also be infiuenced by the symmetry of the order 
parameter and also by other factors [10, 12-14]. In the 
second part of this manuscript we show that specific-heat 
data from DTA measurements can, to some extent, also 
give useful information on absolute values of Cp, and 
in particular reveal trends even in the absence of a 
sharp phase transition. We briefly discuss our data on 
the magnetic field dependence of Cp{H) for VsSi, and 
LuNi2B2C at a fixed temperature T = 7.5 K. 



SAMPLE CHARACTERIZATION 

The VsSi single crystalline sphere (mass m « 11.1 mg) 
used for this study has been previously characterized us- 
ing magnetization measurements by Kiipfer et al. (sam- 
ple "SA" in [15]). The authors irradiated VsSi samples 
with fast neutrons, and the present sample was further 
annealed for 2 h at 630 °C in order to decrease pinning 
effects. We determined the transition temperature to 
superconductivity of the crystal in zero magnetic field 
calorimetrically as Tco ~ 16.8 K. The crystal undergoes 
a martensitic phase transition from cubic to tetragonal 
upon cooling through a temperature Tm ~ 16.7 K (see 
Fig. 3 and later Fig. 6) as soon as superconductivity is 
suppressed by a magnetic field. Fig. 1 shows the phase di- 
agram of the investigated sample, showing the upper crit- 
ical field Hc2 (T) and the weakly field dependent marten- 
sitic transition at Hm{T). 

The LuNi2B2C single crystal investigated here (m w 
13.2 mg) was irregularly shaped, with only a few smooth 
surfaces. No previous data of physical properties were 
available for this crystal. We determined the transition 
temperature as T^q ~ 15.5 K. In Fig. 2 we show the phase 
diagram of this sample, with an unusual positive curva- 
ture of Hc2{T) near Tcq. This is an indication for a high 
quality of the crystal in which the clean limit is achieved 
[16]. 

The 2II-NbSe2 sample under study (m « 6.4 mg) 
consisted of 6 thin slices cut from a large single crystal, 
which we stuck together with Apiczon N grease. The 
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FIG. 1: Phase diagram of the investigated VaSi single crystal, 
showing the upper critical field Hc2{T) (filled circles), and the 
martensitic transition Hm{T) (stars). The magnetic field was 
oriented parallel to the [100] direction. 
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FIG. 2: Upper critical field Hc2{T) of the investigated 
LuNi2B2C single crystal. The magnetic field was oriented 
perpendicular to the layers. Note the unusual positive curva- 
ture of Hc2{T) near Tco- 



sample appeared to be somewhat inhomogeneous since 
we observed two superconducting phase transitions in a 
narrow temperature interval around T = 6.8 K, separated 
by « 80 mK (see below). This inhomogeneity might be 
a property of the crystal or was later induced by the 
cutting process. We show results on this ensemble of 
crystals to demonstrate that high-resolution specific-heat 
data can be used to clearly identify different supercon- 
ducting phases even if the respective critical tempera- 
tures are very close to each other. 
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FIG. 3: Reduced specific heat Cp/T of a VaSi singie crystal 
for four different magnetic fields, measured in a commercial 
PPMS (Quantum Design) by a conventional lieat-relaxation 
technique. 
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FIG. 4: Magnified view of Cp/T of VaSi in zero magnetic field 
in the transition region around Tco- DTA data : open circles. 
PPMS data : squares. The approximate transition width 
(10-90% criterion) for the DTA measurement is marked with 
vertical lines. 



RESULTS 

1. Phase transition in VsSi, LuNi2B2C, and 
2H-NbSe2 

a) V-iSi 

In the first part of this section, we will compare our 
zero-field data acquired with the DTA method with data 
collected in a commercial PPMS (Quantum Design) by 
a conventional heat-relaxation technique. Fig. 3 
shows the temperature dependence of the specific heat of 
the VsSi crystal for four different magnetic fields, mea- 
sured with the DTA technique. While the transition to 
the superconducting state is strongly field dependent, 
the martensitic transition shows only a weak field de- 
pendence (see Fig. 1). The martensitic transition at Tm 
takes place very close to the superconducting transition 
Tco in zero field and therefore causes a certain broaden- 
ing of this transition. 

In Fig. 4 we show both DTA and PPMS data in zero 
magnetic field at the transition region around Tco, mea- 
sured on the same VaSi crystal. Applying a common (10- 
90%)-criterion to the PPMS data, we find AT sa 280 mK 
for the width AT of the phase transition. Near Tco, 
the data-point density for this kind of measurement is 
« 10 data points per Kelvin, and the PPMS setup needed 
about 40 min to measure the displayed data. The reason 
for this relative low data-acquisition rate is that both 
the stabilization of the temperature and the subsequent 
relaxation for every single data point take considerable 
time. In the same figure we show corresponding specific- 
heat data collected with the DTA method on the same 
sample. Near Tco, the data-point density is 165 data 
points per Kelvin, and the DTA experiment took about 
8 min to collect these data. Using again the (10-90%)- 



criterion we find AT « 80 mK for the width of the phase 
transition to the superconducting state. This transition 
is much narrower in the DTA data than in the relaxation 
method data from the PPMS (by a factor sa 3.4). Since 
both measurements have been conducted on the same 
VsSi crystal, this discrepancy in the transition width can- 
not be attributed to sample-property issues, but must be 
ascribed to the different measuring techniques. While ev- 
ery single temperature reading in a DTA experiment can 
represent, in principle, one Cp data point [7, 17], it is the 
temperature interval needed to measure the relaxation 
rate that determines the temperature resolution in a re- 
laxation method experiment. This quantity is usually of 
the order of 100 mK or more, and can only be reduced 
at the cost of a larger uncertainty and scatter in the Cp 
data. 

We next want to compare the DTA data to results of 
corresponding specific-heat measurements published by 
other groups for the same compound that have been ob- 
tained by different methods (see Fig. 5). It is obvious 
from Fig. 5 that the difference in the respective transi- 
tion widths is prominent. We have reproduced the data 
of Sebek et al. [18] in Fig. 5a (open squares) that have also 
been measured using a commercial PPMS system. Here, 
the data-point density is « 4 data points per Kelvin, 
roughly 40 times less than in our DTA measurement. 
The transition width appears to be about AT « 520 mK, 
which is about « 6 times larger than the transition width 
measured with our DTA method. At the same time this 
width is also about twice as large as the transition width 
that we measured for our VaSi crystal with the same 
relaxation method, which may indicate different sample 
qualities. One can safely state here that the relaxation 
methods achieve a considerable lower resolution on the 
temperature scale than what we routinely reach with a 
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FIG. 5: Specific heat of VsSi around Tco in zero magnetic 
field, measured on different samples with different techniques, 
(a) DTA method: single crystal, open circles (this work). Re- 
laxation method: commercial PPMS, same crystal as DTA 
data, filled squares (this work); single crystal from [18], open 
squares; polycrystalline sample from [19], crosses. Adia- 
batic method: single crystal from [20], stars, (b) Heat-pulse 
method: single crystal from [21], triangles; single crystal from 
[22], filled circles. 



DTA method; while the data-acquisition rate is much 
higher with the latter method. 

Next, we compare our DTA data to correspond- 
ing data from a standard semi-adiabatic heat-pulse 
technique [14]. As an example we have reproduced the 
data of Ramirez [19], taken on a polycrystaUine VsSi 
sample in zero magnetic field, in Fig. 5a (crosses). This 
sample was reported to have a Tco = 16.5 K and a tran- 
sition width of - 8% [19] (i.e., AT w 1220 mK), which is 
about 15 times larger than the transition width that we 
measured with our DTA method on a different sample. 
Besides effects due to the different measuring techniques, 
this huge discrepancy can be attributed to the polycrys- 
talline nature of the sample measured in Ref. [19] that 
might lead to a certain inhomogeneity and a broadening 
of the transition to superconductivity. The data-point 
density of that experiment near Tco is 4 data points 
per Kelvin. 

A similar heat-pulse method was also used by Junod 
et al., who measured the specific heat with a classical 
heat-pulse calorimeter on a VsSi single crystal (m w 



3.6 g) [22]. The transition width can be estimated from 
the data reproduced in Fig. 5b to ~ 430 mK, which is 
about five times larger than in our DTA data. The data- 
point density in this experiment is of the order of 5 data 
points per Kelvin. 

As a last example for the heat-pulse method, we 
have plotted corresponding data obtained on a VaSi sin- 
gle crystal by Brock [21]. The transition width, taken 
from the data plotted in Fig. 5b, is « 300 mK, which 
is about four times larger than in our experiment. The 
data-point density near Tcq is of the order of 6 data points 
per Kelvin. 

Finally, we compare our DTA data to an adiabatic 
method as used by Khlopkin [20]. Details about the 
measuring technique arc not given in Ref. [20] and are, 
to our knowledge, not published. We assume that the 
method is either a heat-pulse based technique or a vari- 
ant of a continuous- heating technique [23]. The corre- 
sponding data of a single crystal of VaSi (m =1.5 mg) in 
zero magnetic field are reproduced from Ref. [20] and are 
plotted in Fig. 5a (stars). The data-point density here 
is w 14 data points per Kelvin, with a transition width 
AT w 350 mK. 

The A15 superconductor V^Si usually undergoes a 
martensitic phase transition. This crystallographic phase 
transition is known to be very sensitive to point-like 
or off-stoichiometric defects. Only high-quality samples 
with a resistivity ratio p = i?(300 K)/i?(18 K)> 25 are 
expected to undergo such a transition that can be re- 
vealed by specific-heat measurements. A Peierls instabil- 
ity merely opens a gap on a portion of the Fermi surface, 
and the associated specific heat should behave qualita- 
tively similar to that of a superconductor [24] . The super- 
conducting gap competes with the Peierls gap, however, 
and therefore superconductivity suppresses the marten- 
sitic transition as can be clearly seen from the absence of 
a feature coming from the martensitic transformation in 
our zero-field data shown in Fig. 4. While this transition 
usually takes place well above Tco in most VaSi crystals, 
the martensitic transition in the here investigated crystal 
takes place slightly below Tcq. 

Fig. 6 shows our specific-heat data of VaSi for both 
the PPMS relaxation method and the corresponding 
DTA data in a magnetic field hqH = ST. The transi- 
tion width of the structural phase transition is in both 
cases AT w 790 mK. This coincidence is expected in this 
case since instrumental broadening is most prominent for 
very sharp features, while the specific-heat discontinuity 
of the martensitic transition is too broad to be affected 
by such effects. 



b) LuNiiBiC 

Fig. 7 shows the DTA data of a single crystal of 
LuNi2B2C [m sa 13.2 mg) in the transition region around 



5 



110 




FIG. 6: Magnification of tfie martensitic transition region 
around Tm of VsSi. Filled circles: specific-heat data obtained 
with a commercial PPMS relaxation method. Open circles: 
DTA data on the same crystal measured in a magnetic field 
fioH = 8T. 



Tco for zero magnetic field. The transition width is 
AT sa 220 mK. The data-point density in that temper- 
ature region is « 130 data points per Kelvin. The in- 
set shows the complete data set. Note the asymmetric 
rounding of the discontinuity at the transition to super- 
conductivity at its low-temperature side, which some- 
what complicates the determination of the transition 
width. At present we do not know whether this rounding 
is an intrinsic property of LuNi2B2C or if it is unique for 
our sample. Since this rounding is nevertheless a rather 
narrow feature (« 200 mK wide), it might be smeared 
out in conventional measurements with lower resolution. 
However, if one does not attribute this feature to the 
transition to the superconducting state, the transition 
becomes much sharper, about « 90 mK. 

In Fig. 8 we have plotted the above discussed DTA 
data together with measurements on LuNi2B2C from the 
literature, which we will discuss in the following. As in 
the case of VsSi, the difference in the respective transi- 
tion widths is prominent. 

From the data of Nohara et al. [25] (open squares in 
Fig. 8), taken on a polycrystalline LuNi2B2C sample us- 
ing a relaxation method [26] we obtain, again using 
a (10-90%)-criterion, a transition width AT « 840 mK, 
around nine times larger than we have measured on a 
single crystal using the DTA method. Again, this may 
be related to the polycrystalline nature of the sample in- 
vestigated in Ref. [25]. Here, the data-point density is 
~ 10 data points per Kelvin near Tcq. 

Corresponding data from another polycrystalline sam- 
ple but measured using a heat-pulse method [27] are 
also plotted in Fig. 8 (filled squares). The transition 
width, AT « 420 mK, is approximately 4-5 times larger 
than in our single-crystal DTA-data, which may again 
be related to sample-homogeneity issues. The data-point 
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FIG. 7: Reduced specific heat Cp/T in the transition region 
around Tco of a LuNi2B2C single crystal in zero magnetic 
field obtained with the DTA technique. The approximate 
transition width is marked with vertical lines. The inset shows 
the complete set of data. 



density in this work is w 12 data points per Kelvin. 

As a last example for LuNi2B2C, we refer to the data 
of Kim et al. [28], who investigated the specific heat of 
a single crystal of LuNi2B2C with an unknown method. 
From their data, reproduced in Fig. 8 with stars, we de- 
duce AT « 410 mK for the transition width, again 4-5 
times larger than the transition width measured with a 
DTA method. The data-point density in that work is w 8 
data points per Kelvin. 



c) 2H-NbSe2 

In this section, we will show that the high resolution 
achieved with the DTA method can be used to reveal 
possible sample inhomogeneities that would otherwise re- 
main unnoticed. In Fig. 9a we reproduced the data of 
Sanchez et al. [29] (squares), obtained on a single crystal 
of 2H-NbSe2 with a relELxation method [30]. The data- 
point density is about 15 data points per Kelvin, and the 
transition width AT ss 70 mK. In the same figure, we 
have plotted the data of Nohara et al. [31] obtained with 
a relaxation calorimeter [25, 26] for comparison (open 
circles). The data-point density here is again « 15 data 
points per Kelvin, and AT w 60 mK. 

In Fig. 9b we show corresponding DTA data that we 
have obtained on a NbSe2 sample consisting of several 
single crystalline pieces stacked together with Apiezon N 
grease. The data-point density is about 190 data points 
per Kelvin in the transition region. The DTA data clearly 
show two steps at the transition to superconductivity, 
which indicates the presence of areas with different tran- 
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FIG. 8; Magnified view of the transition region around 
Tco for Cp/T data obtained on different LuNi2B2C sam- 
ples with different techniques in zero magnetic field. DTA 
method: single crystal, filled circles (this work). Relax- 
ation method: polycrystalline sample [25], open squares. 
Heat-pulse method: polycrystalline sample [27], filled 
squares. The data of Kim et al. [28] measured on a single 
crystal with an unknown technique is plotted with stars. 



sition temperatures. The transition width of the whole 
two-step transition is AT k, 80 mK but we expect that 
the transition width of a homogeneous sample with only 
one Tco would be considerably smaller. This two-step 
transition is still rather sharp and can be clearly resolved 
with the DTA method, while it is questionable if such a 
feature could be resolved using a method with larger in- 
strumental broadening and/or lower data-point density. 

To conclude this paragraph, we want to briefly com- 
ment on advantages and some peculiarities of the DTA 
technique as compared to the heat-pulse and relaxation 
techniques. Apart from advantages concerning speed and 
data-point density (see Table I), the resolution of the 
DTA-Cp data on the temperature scale can be adjusted 
after a measurement has been done. The reason is that 
the DTA technique principally measures changes in en- 
tropy S without any significant instrumental broadening 
effects [17]. While in the other two techniques, the reso- 
lution in temperature is fixed by the finite temperature 
intervals used to obtain one single data point in Cp (e.g., 
the pulse height or the temperature range in which a 
relaxation is measured), the resolution in temperature 
of Cp(T) data obtained by the DTA technique is deter- 
mined by numerical parameters that can be adjusted at 
will during data processing after the actual measurement. 
To be more precise, to obtain specific-heat Cp/T data, the 
derivative of S{T) has to be calculated numerically. This 
results in Cp(T) data with an impressive resolution on the 
temperature scale. However, any attempt to increase the 
accuracy of Cp on the temperature scale is at the expense 




FIG. 9: Specific heat of 2II-NbSe2 in the transition region 
around Tco in zero magnetic field. We compare the data from 
Refs. [29, 31] (a) with a corresponding DTA measurement 
on a stack of NbSe2 crystals, (b) The two-step transition to 
superconductivity is clearly visible (see text). The respective 
transition widths are marked with vertical lines. 



of the accuracy in Cp/T as it was discussed in detail in 
Refs. [7] and [17]. The numerical calculation of SS/ST 
from experimental data introduces a certain broadening 
6T in T given by the temperature interval used to cal- 
culate this derivative. Moreover, the scatter Scp in the 
specific-heat data as calculated from such a procedure 
[17] is inevitably influenced by such a procedure. It has 
been shown that the product 5cp6T is a constant that 
is essentially determined by the total heat capacity and 
by the limiting accuracy with which the temperatures 
can be measured [7, 17]. In other words, any attempt to 
increase the accuracy of Cp in T, e.g., by choosing a nar- 
rower interval to calculate 5S/5T ^ leads to an increased 
scatter 5cp. In order to illustrate this issue, we plotted in 
Fig. 10 corresponding data evaluated with two different 
intervals 5T = 10 mK (open circles) and 5T k. 110 mK 
(thick line) to numerically calculate Cp/T from the same 
zero-field S{T) data on VsSi. It can clearly be seen that 
the increase in 5T (here by a factor « 11) leads to a re- 
duction in scatter 5cp, but sharp features like the second 
order phase transition to superconductivity become more 
broad on the temperature scale. A detailed discussion of 
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FIG. 10: Zero-field Cp/T data on VsSi obtained with the 
DTA method and evaluated with two dilferent intervals 5T = 
10 mK (open circles) and ST ~ 110 mK (thick line). 



this issue can be found in Refs. [7, 17]. 

While the present comparison has essentially been 
restricted to the very common heat-pulse and relaxation 
techniques to measure a specific heat, we want to men- 
tion that modern methods of advanced ac calorimctry 
and continuous-heating techniques can reach at least a 
comparable resolution on the temperature scale as DTA 
measurements. As soon as the temperature amplitude in 
an ac experiment is sufficiently small (as, e.g., in the work 
of Lortz et al. [6] using a thermopile- type arrangement), a 
temperature resolution of 1 mK is routinely achieved. As 
far as the data-acquisition time and the data-point den- 
sity is concerned, continuous-heating methods [23] share 
the same advantages with the DTA technique because 
both methods operate with a continuously varying tem- 
perature of the sample platform while monitoring the 
sample temperature as a function of time. 



2. Low-temperature magnetic-field dependence of 
tlie specific lieats of VsSi and LuNiaBaC 

In the last section of this paper we show that specific- 
heat data from DTA measurements can reveal trends for 
the absolute value of Cp even in the absence of a sharp 
phase transition. As an example, we focus on the low- 
temperature specific heats of VaSi and LuNi2B2C that 
smoothly increase with the applied magnetic field H , but 
still may contain valuable information about the physics 
of the mixed state of these compounds. 

In the standard theory of superconductivity [32-34] 
the magnetic-field dependence of the specific heat in the 
mixed state is nearly linear in H for T — > (i.e., pro- 
portional to the number of vortices), and can be written 
as 



Method 


Reference 


Substance 


Single crystal 


Ar(mK) 




DTA 


This work 


VaSi 


Yes 


80 


165 


Relaxation 


This work 


VaSi 


Yes 


280 


10 


Relaxation 


[18] 


VsSi 


Yes 


520 


4 


Heat-pulse 


[19] 


VaSi 


No 


1220 


4 


Heat-pulse 


[22] 


VaSi 


Yes 


430 


5 


Heat-pulse 


[21] 


VaSi 


Yes 


300 


6 


Adiabatic 


[20] 


VaSi 


Yes 


350 


14 


DTA 


This work 


LuNizBaC 


Yes 


220 


130 


Relaxation 


[25] 


LuNiaBaC 


No 


840 


10 


Heat-pulse 


[27] 


LuNiaBaC 


No 


420 


12 


Unknown 


[28] 


LuNiaBaC 


Yes 


410 


8 


DTA 


This work 


NbSea 


6 shces 


< 80 


190 


Relaxation 


[29] 


NbSea 


Yes 


70 


15 


Relaxation 


[31] 


NbSea 


Yes 


60 


15 



TABLE I: Comparison of the (10-90%)-transition widths AT 
of the transition to superconductivity measured with different 
techniques on various samples, pdp is the data-point density 
(in data points per Kelvin). The time needed to take one data 
point is « 3 s with the DTA method, while we estimate it to 
be of the order of minutes for all the other techniques listed 
here. 



Acp/T{H) = {cp{H) - Cp{H = 0))/T = ^nH/Hc2{T) , 

(1) 

where is the Sommcrfcld coefficient characterizing the 
specific heat of the electronic system in the normal state. 
This increase in Acp/T(H) can, in the simplest picture, 
be ascribed to the contribution of the normal cores of the 
flux lines entering the sample, and a low-temperature 
linear term is expected to be proportional to H up to 
Hc2 [35]. At temperatures T > 0, the field dependence 
of Acp/T{H) is expected to be more complicated, and 
a non-linear behavior has been observed in a number 
of compounds: VgSi [19], LuNiaBsC [25, 31], YNiaBaC 
[31, 36], NbSes [29, 31], UPtg [14], CeRus [37], Lu2Fe3Si5 
[38], and in YBaaCugOy [10, 12]. 

We have measured the specific heat of VaSi in mag- 
netic fields up to /io-ff = 8.3 T. Fig. 11a shows two 
DTA data sets for = T and 6 T. Below the tran- 
sition to superconductivity and the martensitic transi- 
tion the two data sets do not match, which is expected 
due to the field dependent Acp/T{H). In Fig. lib we 
show corresponding data in a small temperature interval 
around T = 7.5K for selected magnetic fields. The non- 
linear field dependence at the lowest investigated mag- 
netic fields is already visible in this figure. 

Extracting the Cp/T[H) data for a fixed temperature 
T = 7.5K results in the data plotted in Fig. 12a, demon- 
strating the magnetic-field dependent specific-heat differ- 
ence Acp/T{H) of VsSi (open circles). The dashed line is 
a guide to the eye and visualizes an approximately linear 
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FIG. 11: Cp/T-DTA data of a VsSi single crystal : (a) for 
fioH = OT and 6T; (b) in a small temperature interval 
around T — 7.5 K, starting with ^oH = OT (open circles) 
for the lowest curve and /loH = 4T for the uppermost curve, 
H increasing in intervals of 0.5 T. 



behavior above hqH > 2 T. It can clearly be seen that 
the data deviate from linearity below (j-qH < 2 T where 
they are better described by a iJ^/^ dependence. For 
comparison, we have plotted in the same figure the data 
of Ramirez [19] from a non-transforming polycrystalline 
VaSi sample, i.e., from a sample that does not exhibit 
martensitic transformation. The deviation from linear- 
ity becomes significant for low fields close to the lower 
critical field Hd- In analogy to Rcf. [19] we can fit our 
data to Acp/T{H) = A{H - H*dY/^ with H*^ = 0.2 T 
being the field where magnetic flux starts to enter the 
sample [39, 40], and A w 8.4 mJ mole'^ K'^T"'^/^, in 
good agreement with the result of Rcf. [19]. 

In Fig. 13 we show the magnetic field dependent Cp/T 
data for LuNi2B2C in a similar way as we did for VsSi 
in Fig. 11. From these data we obtain again the specific- 
heat difference Acp/T{H) at a fixed T = 7.5 K. It can 
clearly be seen in Fig. 12b that the data deviate from lin- 
earity below hqH « 1 T where they again may be better 
described by a iJ^/^ dependence. A tentative fit accord- 
ing to Acp/T{H) = A{H - Hci)^/^ with Hci = 87 mT 
[41] gives A « 5.2 mJ mole"! K-^T'i/^. 
For comparison we have reproduced the data of Nohara 
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FIG. 12: Magnetic-field dependence of the specific-heat differ- 
ence Acp/T{H) at r = 7.5 K for different crystals. The thick 
lines represent fits according to Acp/T(H) = A{H — H2\Y^^ 
(see text). The dashed lines are guides to the eye to visu- 
alize the approximate linear behavior for the magnetic fields 
HoH > 2T in the case of VsSi and fjLoH > 1 T in the case of 
LuNi2B2C. (a) Data of a VaSi single crystal obtained from 
DTA measurements (open circles), in comparison with the 
data of Ramirez [19] (filled circles), (b) Data of a LuNi2B2C 
single crystal obtained from DTA measurements (open cir- 
cles), together with the data of Nohara et al. [25] (filled cir- 
cles). 



et al. for LuNi2B2C in Fig. 12b. Our data clearly confirm 
the observation that also in LuNi2B2C Acp/T{H) is not 
proportional to H. 

The prediction of linearity of Acp/T{H) in H strictly 
holds, within the ECS theory, only for T 0. To test 
whether or not the BCS prediction for the behavior of 
Acp/T{H) data taken at fixed but elevated tempera- 
tures might also account for the observed ^/H-like de- 
pendence of Acp/T{H) we have simulated correspond- 
ing data according to a model that assumes a mixture 
of a normal-state electronic specific heat Cen that is ex- 
plicitly proportional to H/Hc2, and a superconducting 
component Ces proportional to (1 — H/Hc2)- For the 
normal-state electronic specific heat we used Cen = InT 
and for the superconducting component a standard BCS 
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FIG. 13: Cp/T-DTA data of a LuNi2B2C single crystal: (a) 
for = OT and 2.5 T; (b) in a small temperature interval 
around T — 7.5 K, starting with /io// = OT (open circles) for 
the lowest curve and ^i^ll = 4T for the upper most curve, H 
increasing in intervals of 0.5 T. 



s-wave expression — 9.17jNTcexp{—1.5Tc/T), which 
is a fair approximation not too close to Tc or Hc2, re- 
spectively. Taking our measured values for Tco and 
T,{H) = T,q{1 - H/H,2y^^ with tioH,2 = 21 T for 
VaSi and 9T for LuNi2B2C, we obtain the differences 
Acp/T{H) as plotted in Fig. 14. It can clearly be seen 
that although the data for T = 7.5 K are somewhat 
rounded as H approaches the upper critical field Hc2{T), 
the distinct peculiar curvature as observed in the corre- 
sponding experimental data for iJ — > (sec Fig. 13) can 
in no way be reproduced. A systematic study of Nohara 
et al. [25] on LuNi2B2C has indeed shown that the ^/~H 
dependence of Acp/T{H) for T is qualitatively pre- 
served up to T « 9 K. 

These findings are noteworthy since both VaSi and 
LuNi2B2C are supposed to be s-wave systems [42-44]. 
The deviation of Acp/T{H) from linearity and especially 
its approximate H^^^ dependence is very often inter- 
preted as an indication for lines of nodes in the energy 
gap and taken as a hallmark for a d-wave symmetry of 
the order parameter [10, 12, 13]. We can state here that 
a non-linear Acp/T{H) as one observes in the s-wave su- 
perconductors VsSi and LuNi2B2C is likely a more gen- 
eral phenomenon of type II superconductors. Possible 
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FIG. 14: Normalized differences l\Cp/'yMT{H) as modeled for 
VgSi (upper panel) and LuNi2B2C (lower panel) for different 
temperatures T, assuming a field dependent mixture of nor- 
mal and superconducting components that vary proportional 
to H/Hc2 and 1 — H/Hc2, respectively (see text). Thick solid 
line: data for T = 0. 



plausible explanations for this peculiar behavior may in- 
clude vortex- vortex interactions [19, 33, 45], the field de- 
pendent shrinking of the vortex cores [31, 46], and multi- 
gap scenarios [47, 48]. 

In summary we have presented DTA-specific-heat 
measurements on various superconductors and compared 
the achieved temperature resolution to corresponding 
data obtained by other methods from the literature. The 
DTA method allows for an exceptionally high resolution 
on the temperature scale while providing a very high 
data-point density and a short measuring time. We have 
also investigated the magnetic-field dependence of the 
specific heats of VsSi and LuNi2B2C at a fixed tem- 
perature T = 7.5 K, and found an excellent agreement 
with corresponding data from the literature. These mea- 
surements demonstrate that low-temperature DTA tech- 
niques are not only suitable to detect sharp phase transi- 
tions, but also allow for sufficiently precise measurements 
of absolute values of c^. 
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